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Shock- triggered formation of magnetically-dominated 
clouds. 

"o II. Weak shock-cloud interaction in three dimensions. 



o 

(N 



O 

m 
< 

6 



> 

oo 
in 

cn 
o 
o 



. Van * S. A. E. G. Falle^ and T. W. Hartquist^ 

School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK 
Department of Astronomy, University of Florida, Gainesville, Florida 32611, USA 
Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK 



Accepted - . Received - ; in original form ■ 



ABSTRACT 

To understand the formation of a magnetically dominated molecular cloud from an 
atomic cloud, we study the interaction of a weak, radiative shock with a magnetised 
cloud. The thermally stable warm atomic cloud is initially in static equilibrium with 
the surrounding hot ionised gas. A shock propagating through the hot medium then 
interacts with the cloud. We follow the dynamical evolution of the shocked cloud with 
a time-dependent ideal magnetohydrodynamic code. By performing the simulations in 
3D, we investigate the effect of different magnetic field orientations including parallel, 
perpendicular and oblique to the shock normal. We find that the angle between the 
shock normal and the magnetic field must be small to produce clouds with properties 
similar to observed molecular clouds. 
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1 INTRODUCTION 

Mol ecular clouds exhib it a hierarchical density structure 
(e.g. iBlitz k Star'^ll986l ) . Stars form in the densest regions, 
or dense cores, which become gravitationally unstable. In 
fact, molecular clouds in the Solar neighbourhood that don't 
harbour any stars are rare. While most of the stars within 
molecular clouds are young 1-2 Myr), stellar associa- 
tions older than 5 Myr are devoid of molecular gas (e.g. 
iBallesteros-Paredes k, HartmannI [2007h . This suggests that 
molecular clouds are short-lived, transient objects and that 
the time-lag between cloud formation and stellar birth is 
short. 

The rapid onset of star formation requires that large 
density contrasts arise while the parental cloud forms. 
It has been suggested that this fragmentation results 
from th ermal processes in the interstellar medium (ISM) 
(see e.g. iKlessen. Krumholz. Heitsch|[20Q9l . and references 
therein). For a range of pressures and the heating and 
cooling rates appropriate for diffuse atomic gas, two ther- 
mally stable phases exist, i.e. a rarefied, warm phase and 
a cold, d ense phase, which can co- ex i st in pressure equi- 
libriu m (iField. Goldsmith, Habind Il969l : IWolfire et all 
Atomic gas at intermediate temperatures is subject 
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to a thermal instability. A sufficient rise in pressure causes 
the cold phase to be the only stable one (lFieldlll965l ). 

Previous studies of molecular cloud formation due 
to thermal instability mainly focus on collisions of warm 
gas streams in the context of expanding supe r novae 
shells or spiral arm shocks (e.g. Hennebelle et al. 2 0081 : 
Heitsch, Stone, & Hartmann 2009; Inoue & Inutsuka 20091) . 
Beside being thermally unstable in some circumstances, the 
collision region is prone to numerous dynamical instabili- 
ties such as the Kelvin-Helmholtz, Ray leigh- Taylor and non- 
linear thin shell or Vishniac instabilities. In the turbulent 
shocked layer cold gas clumps then arise on short timescales. 
Although the derived density and velocity structures depend 
strongly on the magnitude and orientation of the magnetic 
field, they resemble t he ones observed in molecular clou ds 
and diffuse HI clouds (|Heitsch, Stone, HartmannI l2009l V 

While this model of flow-driven structure generation 
shows rapid onset of star formation while the parental cloud 
is forming, it cannot account for the observed low star for- 
mation rates (|Klessen, Krumholz, Heitschll200"9l ) . As there 
is a continuous instream of gas into the collision region, 
too much of the accumulated gas will be converted into 
stars. For the same reason, it also cannot explain cloud life- 
times. Some of these limitations do not necessarily occur 
when studying the same processes in cloud-clo ud interac- 
tions (jKlein k Woodsl Il998l : iMiniati et all Il999l ) or shock- 
triggered models where shocks overrun warm, diffuse density 
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pertu rbations (jlnutsuka Kovamal l2004l : I Van Loo et alJ 
I2OO7I , hereafter Paper I). Indeed, numerical models of shocks 
interacting with clouds s how that the clouds fr agjment due to 
dynamical instabilities ([Mac Low et al.lll994l ). thereby set- 
ting a limit on both the cloud lifetime and the star formation 
rate. 

Colliding flow-driven models only form thin filamen- 
tary clouds, while different cloud morphologies such as a 
cometary cloud structure wi th a massive head an d long- 
spread tail, are also observed (|Tachihara et al.ll2002l V Such a 
morphology is characteristic for clouds harbouring cluster- 
forming cores. The W3 Gian t Molecular Cloud (CMC) is 
an example of such a cloud (|Moore et al.ll2007l V As there 
is no constraint on the geometry of density perturbations 
in the shock-triggered models several morphologies can be 
reproduced. 

The interaction of a shock with a magnetised cloud has 
been studied extensively in a two-dimension al (2D) axisym- 
metric gjeometry, both for a diabatic (e.g. iMac Low et al.l 
I994I: Nakamura et al.l I2OO6I ) and radiative shocks (e.g. 
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limited to a single configuration of the magnetic field, i.e. 
the magnetic field is parallel to the symmetry axis and the 
shock normal. To study the effect of the magnetic field orien- 
tation, it is necessary to model the shock-cloud interaction 
in 3D. Althoug;h such simulations h ave been aroun d for some 
time (|Stone Norman 1992; Gr egori et al.ll2000l ). only re- 
cently high enough resolution has been achieved for genera l 
cloud properties to converge (^Shin. Stone, Snvderl 120081 ). 
However, the Shin et al. study focuses on strong, adiabatic 
shocks, while the results of Paper I show that molecular 
clouds most likely form from interactions involving weak, 
radiative shocks. The first 3D simulations of the interaction 
of a radiative shock w ith a magnetised cloud were performed 
bv lLeao et all (|2009l ). although they use a nearly isothermal 
equation of state to simulate strong radiative cooling. 

In this paper we study the interaction of a weak, radia- 
tive shock with a magnetised cloud. We focus on the early 
stages of the evolution before the cloud re-expands and frag- 
ments. In Sect. [2] we describe the numerical code and the 
initial conditions. The results are presented in Sect. (3] while 
we discuss the cloud properties in O Finally, we finish the 
paper with a summary (Sect.[5|). 



2 NUMERICAL MODEL 
2.1 Numerical code 

To solve the ideal magnetohydrodynamics equations we use 
a seco nd-order G odunov scheme with a linear Riemann 
solver (|Fallelll99lh . To ensure that the solenoidal constraint 
is met, a divergence clea ning algorithm is implemented 
in the numerical scheme (|Dedner et al ] |2002l ). The ther- 
mal behaviour of t he diffuse atomic gas as described by 
IWolfire et al.l (|l995l ) is treated through the inclusion of a 
source term in the energy equation. The exact expressions 
for the cooli ng and heating function are given in Paper I 
and ar e from ISanchez-Salcedo. Vazquez-Semadeni. Gazoll 
(|2002l ). As the cooling time scale can be much smaller than 
the dynamical time scale, the equations are stiff. By using an 
exponential time differencing method the numerical scheme 



remains stable for the larger tim e step set by the Courant 
condition (see e.g. lTokmanll2006l ). 



2.2 Initial conditions 

In our models, a quiescent, uniform and spherical cloud of 
radius (RcO 200 pc and number density n = 0.45 cm~^ is 
initially in thermal equilibrium. The gas is in the thermally 
stable warm phase with a thermal pressure p = 2825/c where 
k is the Boltzmann constant. Also, the cloud is in pressure 
equilibrium with the surround ing hot ionised med i um (n 
= 0.01 cm-^). The analysis of iBegelman McKed (Il990l ) 
shows that the hot component of the interstellar medium 
(ISM) is thermally stable due to the subsequent reheating 
by supernovae. Thus, we assume that the surrounding gas 
cools adiabatically. 

While observations show that molecular clouds are 
magnetically dominated with values of the ratio /3 of the 
thermal gas pressure to magnetic pressure of the order 
0.04-0.6, most diffuse clouds hav e weak magnetic fields 
(iCrutcher, Heiles, Trolandl lioosl ). On scales large com- 
pared to the size of clouds, the magnetic and thermal gas 
pressures are actually comparable. Therefore, we adopt a 
uniform magnetic field with a strength such that (5—1. 
(For our models this means that the magnetic field strength 
is roughly 1 /iC.) 

As in Paper I, we consider a steady, planar shock hitting 
the quiescent cloud. The shock is propagating in the negative 
^-direction with the angle between the shock normal and the 
magnetic field either 0^ (parallel model), 15^^, 45^ (oblique 
models) or 90"^ (perpendicular model). The magnetic field 
is in the x — z plane for the oblique models and in the x- 
direction for the perpendicular model. The shock velocity 
through the hot ionised medium Vext is 2.5 times the hot 
gas sound speed a (i.e. the shock sonic Mach number is 2.5). 
As the fast magnetosonic speed c/ lies within the range of 
values a < Cf < V^a, this shock is a weak, fast-mode shock. 



2.3 Computational domain 

The computational domain is —2Rci < x,y < 2Rci, and 
—2.066Rci < z < 3Rci with free- flow boundary conditions 
on all boundaries except on the positive z boundary. There 
we fix the variables to have the postshock flow conditions 
calculated from the adiabatic Rankine-Hugionot relations. 

For models of adiabatic shock-cloud interaction, it is 
necessary to have at least 100 grid points per cloud radius 
Rci to obtain convergence of general cloud properties like the 
shape of the cloud and the rms velocities along ea ch axis 
(|Mac Low et al.1 Il994l : IShin. Stone, Snvderl I2OO8I ). How- 
ever, some cloud properties do change with increasing res- 
olution. These are usually associated with quantities that 
are sensitive to small-scale pr o cesses . For radiative shocks, 
lYirak. Frank, CunninghamI (|2009h show that there is no 
true convergence at 100 cells per cloud radius. Changes in- 
duced by an increasing resolution have a global effect later 
on in the dynamical evolution of the cloud. It is likely that 
convergence is only achieved when the cooling lengths are 
adequately resolved. At the moment, however, it is not fea- 
sible to perform simulations of a shock-cloud interaction in 
3D that attain the required resolution. Therefore, we adopt 



Weak shock-cloud interaction in 3D 3 



a reso lution similar to that used by IShin. Stone, Snvdei] 
(|2008l l for the adiabatic shocks. We use an adaptive mesh 
with 5 levels of refinement with the resolution at the finest 
grid 480^ X 608. Such a resolution corresponds to a physical 
grid spacing of 1.58 pc. As translucent clumps in GMCs have 
length scales of about 3 pc, it is clear that we cannot resolve 
such structures. Our simulations, thus, can only show the 
onset of clump formation. 

As t he cloud is accelerated to move with the post-shock 
flow fe.g. lMac Low et aLllTooi ^. the cloud eventually moves 
off the grid. To avoid this we calculate the density- weighed 
average velocity of the cloud along each different axis at 
every timestep. The density- weighed average of each velocity 
component is defined as 



</>= 



(1) 



where f is Vx, Vy or Vz and C is a scalar which is 1 for cloud 
material and for ambient gas. p, V and Md are the mass 
density, volume, and cloud mass. By performing a Galilean 
transformation of the computational domain, the cloud re- 
mains in the centre of the grid. Furthermore, we are able 
to study the acceleration of the cloud. Note that density- 
weighed averages can also be calculated for other flow vari- 
ables. 



3 DYNAMICAL EVOLUTION 
3.1 Parallel shock 

The shock-cloud interaction for a parallel shock in 3D is 
nearly identical to the evolution studied in the 2D axisym- 
metric simulations of Paper L It is useful to describe the 
dynamical evolution here again, not only to show the dif- 
ferences compared to the 2D results, but also because some 
of the dynamical characteristics are relevant for the perpen- 
dicular and oblique models. 

Figure [T] shows results for the parallel case at two times. 
As the intercloud shock sweeps around the cloud (in a time 
of tcp = 2Rci/vext), a shock is transmitted into the cloud 
and a bow-shock forms in front of the cloud. The transmit- 
ted fast-mode shock has a lower propagation speed than the 
intercloud shock, i.e. Vint = Vext/x^^'^ where x is the den- 
sity ratio of cloud/intercloud gas. Therefore, a velocity shear 
layer forms at the cloud boundary. Because of this slip sur- 
face, a vortex ring develops and sweeps cloud material away 
from the cloud (see Fig. [T]). 

As the fast-mode shock moves through the cloud, it 
compresses and heats the initially thermally-stable warm 
gas so that it ends up being thermally unstable. Figure. [2] 
shows the distribution of mass in the p — n phase space 
at two different times. The time scale for radiative cool- 
ing is shorter than the propagation time of the fast-mode 
shock through the cloud, given by the cloud-crushing time 
tec = Rci/vint- Thus, the gas loses a significant frac- 
tion of its internal energy during the compression. Fur- 
thermore, the magnetic pressure increases behind the fast- 
mode shock. Hence, the value of /S drops below unity in- 
side the cloud. This provides the ideal conditions for the 
generation, by MHD waves, of dense clumps and cores 
(iFalle k Hartauistll2002l : IVan Loo. Falle. Hart a uisti [20061 : 
IVan Loo et aDl2008l V Cold, dense clumps can also form due 



to small per turbations along the u nstable part of the equilib- 
rium curve (jlnutsuka Kovama| [2007) . In our simulation we 
find a small fraction of cloud material on the thermally un- 
stable part of the equilibrium curve after ^ 6 Myr (or O.Ttcc). 
Typical timescales for both processes are a few Myr. As the 
cloud flattens and fragments in about 1.5-2 tec (see Paper I), 
there is ample time for these processes to work. However, the 
numerical resolution of our simulations is insufficient to es- 
tablish whether these processes are the dominant formation 
mechanisms of dense clumps and cores within clouds. As we 
cannot follow this formation process, we stop the simulation 
shortly after tec- This timescale was chosen because the re- 
expansion phase of the cloud then starts. Also, self- gravity 
which is not included in these simulations becomes globally 
important around this time and will affect the subsequent 
dynamical evolution. 

The fast-mode shock is not the only shock to be trans- 
mitted into the cloud. A slow-mode shock is trailing the 
fast-mode shock. As it moves much more slowly than the 
fast-mode shock, it remains close to the boundary of the 
cloud. Behind the slow-mode shock, the magnetic pressure 
decreases and there is nothing that prevents the gas from 
compressing as it cools. Figure [2] clearly shows the rapid 
condensation due to the slow-mode shock. In a few Myr 
the gas behind the slow-mode shock cools to the thermally- 
stable cold phase. For the gas behind the fast-mode shock 
that is not processed by the slow-mode shock, this process 
occurs much more slowly. 

Thus, a dense, cold layer quickly forms at the cloud 
boundary with the highest densities on the upstream parts 
of the cloud where the shock first hits the cloud (see Fig.[T]). 
While the cooling behind the slow-mode shock is thermally 
unstable, the dense shell is subject to Kelvin-Helmholtz and 
Ray leigh- Taylor instabilities, even though the dynamical in- 
stabi lities are mostly suppr essed by the strong magnetic field 
fe.g. IChandrasekharlll96ll V Hence, the shell breaks up into 
dense fragments. The densest clumps have number densities 
of ~ 10"^ cm~"^ and are a few parsec in size. Note that we 
found roughly the same values for the higher resolution 2D 
simulation in Paper I. For lower resolution 2D simulations 
the densities are not so high. However, the clumps of the 2D 
simulations are essentially axisymmetric rings that break up 
into smaller, higher density clumps in the 3D simulations. 
These cold, dense clumps contain several hundreds of so- 
lar masses each and are gravitationally unstable. Thus, such 
clumps are likely precursors of massive stars. 

3.2 Perpendicular shock 

The dynamical evolution of cloud interacting with a perpen- 
dicular shock is in many ways similar to that of one inter- 
acting with a parallel shock. As the intercloud shock sweeps 
around the cloud, a transmitted fast-mode shock propagates 
through the cloud making the cloud material thermally un- 
stable. Figure. [3] shows the distribution of mass in the p — n 
phase space at two different times. The compression by a 
weak perpendicular shock is somewhat smaller than that by 
a parallel shock. The compression ratio changes from 2.70 for 
a parallel shock to 1.87 for a perpendicular shock. The gas 
behind the transmitted shock is therefore only marginally 
cooler and less dense. 

However, contrary to what occurs in the parallel case, a 



4 S. Van Loo et al. 





Figure 2. Distribution of the mass fraction in a phase diagram for the parallel shock model. The left panel shows the phase diagram 
at 2.5 Myr (or just after tcp) and the right panel at 9 Myr (or just after tec)- The values for the mass fraction range between 10~^ and 
10~^. The black line represents the equilibrium curve. 



slow-mode shock does not arise in the upstream parts of the 
cloud. This can be easily seen in Figs. |4]and [5] which show the 
logarithm of the density for the perpendicular case as a func- 
tion of position for different slices and different times. Slow- 
mode waves, and therefore slow-mode shocks, do not prop- 
agate perpendicular to the magnetic field. However, as the 
intercloud shock sweeps across the cloud, a slow- mode shock 
does arise at the sides of the cloud. This slow-mode shock 
is not as strong as in the parallel case. Thus, it does not 



trigger a condensation near the boundary. Figure [3] shows 
that the transition from warm, rarefied gas to cold, dense 
gas is rather smooth with most of the gas near equilibrium. 
As discussed in Sect. 13.11 such a situation can initiate the 
for mation of dense, cold clump s embedded in warm, rarefied 
gas (jlnutsuka k Kovamall2007l V 

The boundary layer in the perpendicular shock model is 
thus not as dense as in the parallel one. The maximum num- 
ber density is an order of magnitude smaller ^ 180 cm~^. 
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Figure 3. Similar to Fig. [2] but for the perpendicular shock model. 



Furthermore, the boundary layer is not subject to strong 
dynamical instabilities. While the boundary fragments into 
small, high-density clumps in the parallel model, it does not 
for the perpendicular model. Dynamical instabilities, no- 
tably the Kelvin-Helmholtz instability, take longer to de- 
velop due to a lower velocity shear and lower density, and 
are suppressed due to the strong magnetic field. 

As seen in Figs. [H and \5\ which are of the logarithmic 
density for different slices, instabilities are generated after 
some time, but only in a specific plane, i.e. in the x = 
plane. It is the only direction for which a shear layer forms 
in the perpendicular model. In contrast, in the parallel shock 
model the shear layer covers the entire cloud. The reason for 
this directional preference is found through the considera- 
tion of the vorticity equation, i.e. 

^ = V X (v X cj) + ^ Vp X V(p + ^) - ^ Vp X [(B.V)B] 

where cj = V x v is the vorticity. While the first term on 
the right hand side describes the stretching of vorticity, the 
second is the baroclinic term and accounts for the changes 
in the vorticity due to the intersection of density and ther- 
mal and magnetic pressure surfaces. The third term relates 
to the changes due to the interaction of the density gradi- 
ent and magnetic tension forces. For flows with /3 ^ 1, the 
generation of vorticity is dominated by the thermal pres- 
sure contribution to the baroclinic term. This term does not 
introduce any directional preference. For the parallel shock 
model, the post-shock flow has a large /3 and vorticity is 
generated all around the cloud. Figure [6] shows the down- 
stream plasma beta as a function of the angle between the 
shock normal and the magnetic field. For flows with /3 ^ 1, 
the contribution of the magnetic field to the vorticity gen- 
eration cannot be neglected. Magnetic tension forces then 
annul the generation of vorticity by the baroclinic term. It 
is thus clear that in regions where the magnetic field lines 
are being bent, no shear layer arises. For the perpendicular 
shock model, this is exactly what happens. The post-shock 



flow behind the intercloud shock drags the magnetic field 
with it along the z-direction. The field lines anchored in the 
cloud are bent and stretched in the x, ^-plane near the cloud 
surface. Hence, less vorticity is generated in the x, 2;-plane 
and a shear layer is only present in the x = plane. 

3.3 Oblique shock 

The perpendicular and parallel shock models describe the 
extreme cases of the shock-cloud interaction. Therefore, we 
can expect that the dynamical evolution for oblique shocks 
lies in between those of the perpendicular and parallel mod- 
els. 

For an angle of 45^ between the shock normal and 
the magnetic field, the evolution of the cloud is similar to 
that of the perpendicular shock model. Although the tran- 
sition of the gas from the warm to the cold phase is not 
as smooth as for the perpendicular shock model, it still pro- 
gresses along the equilibrium curve. The maximum densities 
in the cloud are somewhat higher than for the perpendicular 
shock model. As for that model, the post-shock /3 is ~ 1 (see 
Fig. [6]). Thus, magnetic tension again prevents the formation 
of a shear layer surrounding the cloud. The preferential di- 
rection of shear layer is perpendicular to the initial magnetic 
field and the flow direction. 

The 15^ oblique shock model does not resemble the ge- 
ometrical evolution of the perpendicular shock-cloud inter- 
action. Instead, the evolution better matches that of the 
parallel shock model as expected from the value of /3 of the 
gas downstream the intercloud shock (Fig. [6]). As the inter- 
cloud shock sweeps around the cloud, a shear layer forms 
all around the cloud. However, a higher level of shear is 
found in the plane perpendicular to the magnetic field and 
the flow direction. As for the parallel shock model, there is a 
rapid condensation behind the slow-mode shock and a dense 
boundary layer forms. This dense boundary layer, however, 
does not fragment due to instabilities as the small transverse 
component of the magnetic field has a stabilising effect. Al- 
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Figure 4. Similar to Fig. [T] but for the perpendicular shock model. 





Figure 5. Similar to Fig. 0] but the slice is now in the x = ^ plane. 



though the geometrical features of the cloud resemble those 
of the parallel shock model, the cloud properties resemble 
those of the perpendicular model (see Sect.[3|). 

The shock-cloud interaction can thus be separated into 
two evolutionary tracks, i.e. one that is quasi-parallel and 
one that is quasi-perpendicular. From Fig.[6l we see that the 
transition is associated with the value of (5 downstream of 
the intercloud shock. For parallel shocks the post shock flow 
has a /3 > 1. Hence, the magnetic field is dynamically unim- 
portant in the external gas which allows a vorticity layer to 
form around the cloud. For angles > 20^, the magnetic field 



is important and suppresses the generation of vorticity and 
shear layers. 



4 DISCUSSION OF CLOUD PROPERTIES 
4.1 Size and density 

In order to study the properties of the different shock- 
cloud interaction models, we use diagno stics similar to 
those used by previous author s, (e.g. iMac Low et al.l 
Il994l : IShin. Stone, k Snvderllioosi ). Specifically, we use the 
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Figure 6. Plasma beta of the gas downstream the intercloud 
shock as a function of the angle between the magnetic field and 
shock normal. 



Figure 7. Mean cloud volume derived from the mass-weighted 
moments along the different axes for the parallel (solid), per- 
pendicular (dotted) and oblique shock (thick dashed for 45° and 
dashed for 15°) cases. The volume is normalised to the initial 
cloud size. 



density-weighted average of variables such as the plasma /S 
and the density using Eq. [T] 
Furthermore, we define 

a= [5(< x^>-<x>^)Y^\ 

along the x-axis, with similar expressions for analogous 
quantities defined along the y and z-axes. This gives us the 
axes for an ellipsoid with a similar mass distribution as the 
cloud's which are used to follow geometrical changes of the 
cloud. Figure [71 shows the temporal evolution of the volume 
of the cloud. The parallel shock compresses the cloud the 
most, while the perpendicular shock is the least compressive. 
At tec the volume of the cloud is only a tenth of its original 
volume. The ram pressure exerted by the post-shock flow of 
the perpendicular shock is smaller than that for the paral- 
lel shock, as both the post-shock density and gas velocity 
(in the lab frame) are lower. (A parallel shock compresses 
gas more than a perpendicular shock with the same Mach 
number.) Furthermore, work needs to be done to compress 
the transverse magnetic field. Similarly, it is not surprising 
to find that the 45° oblique shock squeezes the cloud faster 
than the perpendicular shock, but slower than the parallel 
shock. For smaller angles we expect the confinement of the 
cloud to converge to the parallel shock model. Indeed, the 
volume of the cloud in the 15° oblique shock model changes 
as for the parallel shock model with only a small deviation 
at later times. 

Figure [7] seems to show a minimum in the volume near 
9 Myr. From the 2D axisymmetric simulations, we know 
that a minimum arises as the fast-mode shock is reflected 
at the centre of the cloud at around tcc(~ 8.3 Myr) and 
propagates back to the boundary. This shock reflection ini- 
tiates a re-expansion phase in the direction perpendicul ar to 
the post-shock flow direction fsee lMac Low et aLlllQQ^ . Pa- 
per I) Along the flow direction, the compression of the cloud 
continues and the cloud ends up as a thin disc after 1.5-2 
tec- This means that the geometrical shape of the cloud re- 
sembles more closely an oblate spheroid than a sphere. The 
ellipticity of the cloud changes from initially (i.e. a sphere) 
to roughly 0.5 at tec- As mentioned earlier, we do not follow 




time (Myr) 

Figure 8. rms number density of the cloud for the parallel (solid), 
perpendicular (dotted) and oblique shock (thick dashed for 45° 
and thin dashed for 15°) cases. 



this expansion phase and stop the simulation shortly after 
the cloud-crushing Lime tec- 
Figure [8] gives the evolution of the root-mean-square 
(rms) number density. Although the volumes of the cloud 
are only different by a factor of 2 as can be seen from Fig. [3 
the rms density changes more significantly. At tec, the rms 
number density for the parallel model is ~ 40 cm""^, while 
it is only 10 cm""^ for the perpendicular model. This large 
variation is due to the large difference in the density of the 
boundary layer between the models (see Sect. [3]). Only mod- 
els with small angles between the shock normal and the mag- 
netic field have mean densities similar to those of observed 
GMCs. Large angle models only reach mean densities sim- 
ilar to diffuse HI clouds. Such a dependence of the density 
str ucture on the magnetic field orienta tion is also observed 
bv lHeitsch. Stone, HartmannI (l2009l ). 
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Figure 9. Mass- weighted plasma beta of the cloud for the parallel 
(solid), perpendicular (dotted) and oblique shock model (thick 
dashed for 45^ and thin dashed for 15^). 
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Figure 10. Mass fraction of the cloud with /3 < 0.1 for the paral- 
lel (solid), perpendicular (dotted) and oblique shock model (thick 
dashed for 45^ and thin dashed for 15^). 



4.2 Thermal and magnetic pressure 

Observations show that molecular clouds are magneti- 
cally dominated with plasma o f the order 0.04 - 0.6 
(jCrutcher. Heiles. Trolandl l2003l ). Figure [9] shows the 
mass- weighed mean of /3. This weighted mean value of /3 for 
the parallel shock model does not lie within the observed 
range. During the early stages of the evolution the gas be- 
hind the slow-mode shock has a large thermal pressure, while 
the magnetic pressure is small. Hence, the plasma f5 is large 
for the high-density gas in the boundary layer (which domi- 
nates the mean value of /3), i.e. /3 5. After about 2.5 Myr, 
the weighted mean value of f3 decreases as the thermally un- 
stable gas behind the fast-mode shock becomes magnetically 
dominated. Although the weighted mean value of f5 is now 
around unity. Fig. [10] shows that a significant mass fraction 
of the cloud has /3 < 0.1. At tec, about 75% of the gas is 
significantly magnetically dominated. 

The other models do produce a cloud with a weighted 
mean 13 of the order of 0.4 (see Fig. [H). The mam reason 
for this has to do with the role played by the transverse 



component of the magnetic field. While the thermal pres- 
sures behind the fast-mode and slow-mode shock are not as 
high as in the parallel shock model (i.e. the transition from 
the warm phase to the cold phase is smoother) , the magnetic 
field is significantly compressed behind the fast-mode shock. 
The combined result of these effects produces a lower (5 in 
both the boundary layer and within the cloud. Figure [10] in- 
deed shows that magnetically-dominated gas appears much 
earlier for shock models with a transverse component to the 
magnetic field. After only 3 Myr, already 10% of the total 
mass in these models is magnetically dominated. 

Our simulations show that, to produce clouds with 
magnetically-dominated high-density gas, the angle be- 
tween the shock normal and the magnetic field must be 
small. While perpendicular shocks produce magnetically- 
dominated gas at low densities, high-density clumps with 
/3 ^ 1 arise in the parallel shock models. 



4.3 Fragmentation 

Our models show that large fractions of the cloud are 
magnetically dominated. This provides the ideal condi- 
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shock models with a transverse component of the magnetic 
field, this process initiates earlier suggesting a higher degree 
of fragmentation. On top of that, another process is effective 
in these models. As the transition from the thermally-stable 
warm phase to the cold one follows the unstable part of the 
equilibrium curve, small perturbations can initiate the for- 
mation of dense, cold clum ps embedded in warm, diffuse gas 
(jlnutsuka Kovamall2007l l . Unfortunately, we cannot follow 
these clump and core formation processes as our resolution 
is insufficient. 

While a low resolution is partly to blame for the 
low amount of fragmentation, the uniform initial condi- 
tions of the cloud also play an important role . In th e col- 
lidin^ ffow-driven models of Hennebelle et al.l (|2008l l and 
iHeitsch. Stone, h Hartmannl ( 2009l V the generation of cold 
dense cores and clumps relies on seeded perturbations in 
either the incoming fiow or at the collision front. With- 
out these perturbations, the collision region remains roughly 
uniform and fragmentation occurs on very long timescales. 
Therefore, it can be expected that the introduction of per- 
turbations within the cloud and at the edge of the cloud 
would produce much more fragmentation. 

Furthermore, our simulations do not include the effect 
of self-gravity. While self-gravity is dynamically unimpor- 
tant for the global evolution of the cloud, i.e. the mass of the 
cloud is much smaller than its Jeans' mass, its effect will be- 
come important locally. For example, the high-mass bound- 
ary layer clumps of the parallel shock model have masses 
that exceed their Jeans mass. We will investigate the effect 
of self- gravity on the cloud evolution in a later paper. 



4.4 Molecular clouds 

Our simulations only describe the dynamical evolution of a 
cloud from warm atomic gas to cold atomic gas. Figure [TT] 
shows the temporal evolution of the cold gas mass fraction 
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Figure 11. Mass fraction of cold (solid) and molecular (dashed) 
gas within the cloud for the parallel (thick) and perpendicular 
(thin) shock model. Molecular gas is assumed to be present when 
the thermal gas pressure is higher than 2500k and the gas tem- 
perature lower than lOOK. 

for the parallel and perpendicular shock model. Cold gas 
arises earlier in the parallel shock model than in the per- 
pendicular one, but after 9 Myr more than 50% of the ini- 
tial cloud mass is in the thermally stable cold gas phase. 
Although we have not included a description of molecular 
cooling, nor do we follow the cloud chemistry, we can roughly 
estimate how muc h ^as is converted from th e cold atomic gas 
to molecular gas. IWilUams. Blitz, Star^ (|l995l ) find that 
the average number density of H2 in the CO clumps of the 
Rosette Molecular Cloud is ^ 220 cm~^. With typical exci- 
tation temperatures between 10 and 20 K, the thermal gas 
pressure of these clumps is roughly 2500k. Therefore, we as- 
sume that any gas parcel in the shocked cloud with a ther- 
mal pressure higher than 2500k and a temperature below 
lOOK will become molecular given enough time (see below). 
Using this criterion, we find that the parallel shock model 
generates a molecular cloud as half of the cold gas becomes 
molecular (see Fig. [TT)) . In the perpendicular cloud model 
only a small fraction of the cold atomic gas is converted into 
molecular gas. This suggests that the perpendicular shock 
model only produces a diffuse HI cloud. 

The above result is only valid if the time scale for the 
formation of molecules is short and if the physical param - 
eters of the formation process are met. lClover et al.1 (|2009l l 
use high-resolution 3D simulations of turbulent interstellar 
gas to follow the formation and destruction of molecular 
hydrogen and CO. They find that most CO forms within 2- 
3 Myr for dense, turbulent gas, while the formation of H2 is 
even faster, i.e. within 1-2 Myr. Their results indicate that 
once large enough spatial and column densities are reached, 
the conversion from atomic to molecular gas is rapid. A good 
indictor for the formation of molecules i s the visual extinc- 
tion Ay, which can be expressed as fe.g. IChapuis Corbell 
I2OO4I ) 



1.80 X 102icm-2' 

For regions with Ay ^0.5 and high local densities, we can 
then expect that molecules are present. Figure [12] shows that 
the visual extinction is already high early on in the parallel 



shock model. These regions also correspond to high-density 
regions, as can be seen in Fig.[T] This model thus most likely 
produces a molecular cloud. Also, note the similarity of the 
column density plot for the parallel shock model with the 
emission map of the W3 CMC (see Fig. 6 of Paper I). A more 
structured inner cloud can be expected with a non-uniform 
initial condition and a higher resolution. 

While the 15° oblique shock model also produces high 
column densities which coincide with high density regions, 
the perpendicular and the 45° models do not. Models with 
large transverse components of the magnetic field produce 
diffuse HI clouds instead of molecular clouds. However, this 
conclusion only holds for our current simulations. A higher 
resolution and inclusion of small-scale perturbations poten- 
tially would produce higher density clumps for these models. 



5 SUMMARY 

In this paper we have presented 3D simulations of the inter- 
action of a weak, radiative shock with a magnetised, diffuse 
atomic cloud. The interaction of the shock induces the tran- 
sition of the cloud from the thermally warm atomic phase 
to the cold one. By modelling the shock-cloud interaction in 
3D, we are able to study the effect of different magnetic field 
orientations including parallel, oblique and perpendicular to 
the shock normal. 

Contrary to the strong, adiabatic shock models of 
IShin. Stone, Snvderl (I2OO8I I, we find that the structure 
of the shocked cloud differs significantly with the magnetic 
field orientation. The shock-cloud interaction can be sepa- 
rated into two distinct classes, i.e. a quasi-parallel one and 
a quasi-perpendicular one. In the quasi-parallel shock mod- 
els high-density clumps are generated in the boundary layer 
surrounding the cloud. As the visual extinction of the gas 
is also high, the resulting cloud is most likely molecular. 
The quasi-perpendicular shock models, however, only pro- 
duce low-density clouds resembling HI clo uds. Thi s resul t is 
similar to the one of Heitsch, Stone, & H artmannI (|2009l l. 

All our models show that the shocked cloud be- 
comes magnetically dominated after a few Myr. Al- 
though this provides the ideal c onditions for the forma- 
tion of dense clumps and cores (iFalle Har tquist 20o3: 
IVan Loo. Falle. Hartauistll2006l : IVan Loo et^a l. 2008 ), we 
do not see this happening in our simulations. This can 
be partly ascribed due to the assumption of an initially 
quiescent, uniform and spherical cloud. From the colliding 
flow-driven models we know that the inclusion of small- 
scale perturbations will generate a higher degree of struc- 
ture. Increasing the resolution can have a similar effect 
(jYirak, Frank, Cunninghamll 20091 ) . In a subsequent paper 
we will investigate the effect of small-scale perturbations and 
a higher resolution on the formation of clumps and cores. 
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Figure 12. Visual extinction along the y-axis for a parallel shock interacting with a diffuse cloud. The left panel show the visual 
extinction at 2.5 Myr and the right panel at 9 Myr. 
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